Relativistic polytropic spheres embedded in a chameleon scalar field 



Vladimir Folomeev, 1,2 * Douglas Singleton 3 ' 

1 Institut fur Physik, Universitdt Oldenburg, 
Postfach 2503 D-26111 Oldenburg, Germany 
2 Institute of Physicotechnical Problems and Material Science of the NAS of the Kyrgyz Republic, 
265 a, Chui Street, Bishkek, 720071, Kyrgyz Republic 
3 Physics Department, California State University, 
Fresno, Fresno, California 93740-8031, USA 



In this paper we investigate gravitationally bound, spherically symmetric equilibrium configura- 
tions consisting of ordinary (polytropic) matter nonminimally coupled to an external chameleon 
scalar field. We show that this system has static, regular, asymptotically flat general relativistic 
solutions. The properties of these spherical configurations, such as total mass, distribution of mat- 
ter, and size, depend strongly on the surrounding scalar field. The mass is found in terms of the 
parameter a, which is proportional to the central mass density of the ordinary matter. We perform 
a stability analysis of these spherical solutions and find an upper bound for a where dynamical 
instability occurs. 

PACS numbers: 04.40. Dg, 04.40.-b, 97.10.Cv 

I. INTRODUCTION 

Scalar fields have been applied to a wide range of model building in cosmology and astrophysics. At the cosmological 
scale scalar fields have been used to model inflation and the late time acceleration due to dark energy. At the galactic 
scale scalar fields have been tried as models for dark matter (for a review with extensive references to scalar field 
models of dark energy and dark matter, see [l|). At the astrophysical scale there are boson stars (see the reviews 
Boson stars are stellar size or smaller, self-gravitating configurations of scalar fields. The scalar fields may 
be self- interacting or not, and one can include a fluid which interacts with the scalar field either only gravitationally 
[3, 0] or through direct coupling Q . The properties of these boson stars strongly depend on the type of scalar fields 
of which they are composed as well as how these fields interact with each other and the surrounding matter/fluid. 

In scalar-tensor gravitational theories there is also the possibility of constructing mixed configurations of a scalar 
field plus ordinary matter. In such theories, a conformal transformation from the original Jordan frame to the Einstein 
frame is always accompanied by the appearance of a nonminimal coupling between ordinary matter and a scalar field, 
which leads to new features not present in general relativity. In particular, this nonminimal coupling between the 
matter of a neutron star and the scalar field can lead to the effect of spontaneous scalarization of the neutron star 
In turn, this leads to a redistribution of the scalar field inside the neutron star, which has a significant influence 
on the process of the neutron star's collapse @. 

In cosmology, this ability of scalar fields to interact nongravitationally with other fields has been used within the 
framework of chameleon cosmology [Toj - TT^ ] , where a chameleon scalar field (j> interacts directly with ordinary matter 
through a conformal coupling of the form e°"^. The name "chameleon scalar field" was suggested in [l(| since the 
properties of the scalar field, such as its mass, depend sensitively on the environment in which the scalar field finds 
itself. The idea that the properties of a scalar field could be influenced by the environment/matter surrounding the 
scalar field was studied earlier, in particular, in the papers [l3l.ll4j where the interaction between matter and the scalar 
field was used to model a dependence of fundamental coupling constants on the local environment. The possibility 
of having scalar fields interact directly with ordinary matter has been used in [12| to describe the present accelerated 
expansion of the Universe. These cosmological models were developed further in where the authors choose a 
generalized expression for the nonminimal coupling between the scalar field and the matter having the form f(<fi)L m , 
with /(</>) being an arbitrary coupling function and L m being the Lagrangian of ordinary matter. In this case, for 
certain parameters of the model, the numerical calculations of [151 ] showed that it was possible to describe the evolving 
Universe with a contraction phase, a subsequent bounce phase, and finally a phase of accelerated expansion. One can 
also find analytical solutions [l6| within the framework of the cosmological model of . These analytical solutions 
describe the evolution in time of the chameleon scalar field. 
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It is alsopossible to consider the effect of a chameleon scalar field at astrophysical scales rather than cosmological 
scales. In [17j the influence of a chameleon field on the structure of a star composed of a polytropic fluid nonminimally 
coupled to a chameleon scalar field was studied. It was found that for special choices of the scalar field potential energy 
and the coupling function /(</>), one could obtain static, regular, asymptotically flat solutions in both relativistic and 
nonrelativistic limits. However, a preliminary stability analysis of the solutions found i n II 711 indicated that they were 
not stable in that the matter of the solutions would either collapse or disperse. In [18[ it was shown that taking 
an isothermal fluid coupled to a chameleon scalar field led to both stable and unstable static solutions, and for the 
unstable solutions, the problem of the gravitational collapse of the nonrelativistic configurations has been considered. 

In this paper we continue the study of the influence of a chameleon scalar field on the structure of polytropic spheres 
begun in [17J . We consider a relativistic, polytropic, spherically symmetric fluid embedded in an external chameleon 
scalar field, which is distributed homogeneously and isotropically over the Universe. When polytropic matter is placed 
in such a homogeneous background, it begins to interact with the background scalar field both gravitationally and 
through the nonminimal coupling. This results in a change of the inner structure of the polytropic star depending 
on the properties of the surrounding chameleon scalar field. In general, the scalar field will evolve over cosmological 
time, and this will lead to a time dependence for the term /(0)L m , which in turn can have a strong influence on the 
structure of the compact fluid configuration which is embedded in the scalar field. In this paper we restrict ourselves 
to only a simplified version of the problem, when we just study changes in the structure of a static polytropic, fluid 
sphere embedded in an external dynamic chameleon scalar field. For the time variation of the scalar field we take the 
time dependence suggested in [l6[ (see below in Sec. IIIIBj) . 

The paper is organized as follows: In Sec.|TT]the general equations describing an equilibrium configuration consisting 
of a real scalar field coupled to a perfect polytropic fluid are derived. In Sec. IIIII these equations are written for a 
particular case when the coupling function /(</>) is chosen to be the exponential one. In Sec. MI 51 we give the results 
of the numerical calculations for this choice of nonminimal coupling function. In Sec. IIVI we address the issue of 
the stability of the obtained equilibrium configurations against small radial perturbations. Finally, in Sec. |V] we 
summarize the main results. 



II. DERIVATION OF EQUATIONS FOR HYDROSTATIC EQUILIBRIUM CONFIGURATIONS 



As discussed in the Introduction we consider a gravitating system of a real scalar field coupled to a perfect fluid. 
The Lagrangian for this system is 

c 4 

L = —r-^zR + -dnvd^tp - V(v) + f(v)L m . (1) 



-f-R + \d»<pd"<p - V(<p) + f(fp)L T 



Here <p is the real scalar field with the potential V(<p); L m is the Lagrangian of the perfect isotropic fluid i.e. a fluid 
with only one radial pressure; f(<p) is some function describing the nonminimal interaction between the fluid and the 
scalar field. The case / = 1 corresponds to no direct coupling between the fluid and scalar field, but even in this case 
the two sources are still coupled via gravity. 

The Lagrangian for an isentropic perfect fluid has the form L m = p 
ing energy-momentum tensor is (details are given in the appendix of 



19, 20]. Using this Lagrangian, the correspond- 

3) 



rpk 



f [(e + P ) Ul u k - 5<<p] + di<pd k <p - Si 



d^ip - v((p) 



(2) 



where e and p are the total energy density and the pressure of the fluid, and u l is the four-velocity. Bearing in mind 
that we will consider below small motions of the matter in the radial direction, we take the time-dependent metric of 
the form 



ds 2 



e u d{x ) 2 -e x dr 2 -r 2 dQ 2 , (3) 

where v and A are functions of the radial coordinate r and time coordinate x a = ct, and dfl 2 is the metric on the unit 
two-sphere. The (§), ({) and (J) components of the Einstein equations for the metric §5§ and the energy-momentum 
tensor ((2|) are then given by 



G° = -e 



-x 



1 



G\ 



—A 
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87rG >0 
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8ttG , 








8?rG i 






C 4 V 







(4) 
(5) 
(6) 
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In the above equations, the prime and dot denote differentiation with respect to r and a; , respectively. 
The equation for the scalar field p coming from the Lagrangian ([!]) is 



_L d_ 



-99 



dp 
dx k 



_ _dV L d£_ 

dip m d<p 



(7) 



Not all of the Einstein field equations are independent because of the law of conservation of energy and momentum, 



T* k = 0. Taking the i = 1 component of this equation gives 



8T? 8T} I, ■ n 1 , < n , , 2 



1 



0. 



(8) 



Thus we have five unknown functions: is, A, ip, s, and p. Keeping in mind that e and p are related by some equation 
of state, there are only four unknown functions. In this section we consider static equilibrium configurations so that 
Eq. ([6]) is trivially satisfied, leaving four equations: the two Einstein equations (j4|) and ([5]), the scalar field equation 
(O, and the hydrostatic equilibrium equation Using the energy-momentum tensor the right-hand sides of 
Eqs. and ([5]) take the form 



J 



1 



X, J2 



A + 2 e ^ 



V(p), 



Tl = -fp--e- x p l2 + V{p). 



(9) 
(10) 



In turn, the scalar field equation ([7]) [with the perfect fluid Lagrangian of the form L m — p, and with the metric 
is 



9 



dV df 

dip ^ dip 



(11) 



Next, taking into account the corresponding components of the energy-momentum tensor from ([2|): 

1 



-e- x p' 2 



V(p) 



and Tq, Tl from ©, (fTUf and using pT|) . allows us to write ([5]), in the static case, as 

dp 1 ^ c?^ 
dr 2 dr 



(12) 



Finally, we need an equation of state for the fluid. We use the following parametric relation between the pressure 
and total energy density of the fluid |21j : 



p = Kp\+ l '« 



PgC 



np. 



(13) 



The constant n is called the polytropic index and is related to the specific heat ratio 7 via n = 1/(7 — 1), p g is the 
rest-mass density of the fluid consisting of a gas of noninteracting particles, and K is an arbitrary constant whose 
value depends on the properties of the fluid under consideration. 
Introducing the new variable 9 defined as in [22j . 



Pg = P 9 cV , 

where p gc is the central density of the fluid, we may rewrite the pressure from Eq. (|13|) in the form 



Using JT3J) and IQ2J in ([l2j leads to 



p = Kp]+ 1/n n+1 . 



sdO , , , „. dv 
2a{n + 1)— = - 1 + a(n+ 1)6} — , 
dr dr 



(14) 
(15) 

(16) 



with a = Kpgi™ /c 2 — p c /(p gc c 2 ) and p c is the pressure of the fluid at the center of the configuration. This equation 
may be integrated to give e v in terms of 9: 



' 1 + <r{n + 1) 
l + a(n + l)e 



(17) 



4 



where e Vc is the value of e v at the center of the configuration where 6 = 1. The integration constant v c , corresponds 
to the value of v at the center of the configuration. It is fixed by requiring e v = 1 at infinity i.e. that the space-time 
be asymptotically flat. 

Thus, the gravitating system of a real scalar field interacting with a perfect fluid is characterized by three unknown 
functions - A, 8, and ip. These three functions are determined by the three equations (|4j) , ([5]) and (jlip . and also by 
the relation (|17|) . We rewrite these equations by introducing a new function M(r), 



p -x = 1 _ 2GM{r) 



M(r) is the mass of the configuration in the range < r < i?, where R is the boundary of the fluid where 
Using this function, Eqs. (j4j and ((9|) yield 



dM _ 4vr 
dr c 2 



2GM(r) 



(18) 
= 0. 

(19) 



To avoid a singularity at the origin, one has to require that M(0) = 0. 

Since we will analyze the system of equations Q , (JS|) and (jll[) numerically, we introduce dimensionless variables 



AttG 



a(n+ l)c 4 



1/2 



tp(r), where A 



A-KGpg C 



a(n + l)c 2 



1/2 



where A has the dimensions of inverse length. With this one can rewrite (j4]) and ([5|) in the form 



(20) 



dv 

^2 1 " £ ^ 

s l + fj(n + l)6»^ 



= c 2 U(i + ^)r 



= e 



-M 



n+l 



1 - 2a(n + l) 



l-2a(n+l)- 



where V = V/(p gc c 2 ) is the dimensionless potential energy of the field. 
Next, using ^51) and (flg]>. we rewrite PT|) as follows: 



d? 



2 cr(n + 1) 

f ~ l + a(n + l) 



l-2o-(n+l)- 



d£ 



V 



+v 



d0 l + a(n + l)6 1 f dv v 



dV 
~d4 



di 



(21) 
(22) 



(23) 



Thus, the static configuration under consideration is described by the three equations (|2"Tj) -(l2U |) . 



III. EQUILIBRIUM CONFIGURATIONS WITH / = ae 



Here we use the equations obtained in the previous section with the specific choice of the coupling function /, 

/ = ae- h \ (24) 



where a, b are arbitrary constants. Such a form for the coupling function was used in [15( to study the evolution of 
the Universe within the framework of chameleon cosmology. In this paper it was shown that choosing different values 
of the parameters a, b gave different forms of the effective equation of state of matter filling the Universe. 

Here we restrict ourselves to the simplest variant when the potential energy is absent i.e. V(4>) = 0. In the absence of 
the nonminimal coupling, this corresponds to the case of a massless scalar field. However, in the case being considered 
here one can see from the right-hand side of (|23|) that the presence of the term with a nonminimal coupling leads to 
an effective mass term of the form a0 n+1 df /dep. Our purpose is to study the influence that such a chameleon scalar 
field can have on the inner structure of polytropic spherical configurations. In this case the total energy density of 
the system Tq from (J5]) takes the form (in units of p gc c 2 ) 



T, 



9 - ne ~H 



(1 



l-2a(n+l)- 



d£ 



(25) 
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Because of the spherical symmetry of the problem, one needs to require that d(j>/d£ = at the center. To ensure the 
positiveness of the energy density over the whole volume of the sphere, it is necessary to choose a > 0. Furthermore, 
from (|2T j) - ([23l) the parameter a can be absorbed by introducing the rescaling y/a£ = £ and y/av — v. Thus we set 
a = 1 in further calculations and drop the bar for simplicity. But we will bear in mind that a can always be restored 
by using the above rescaling, leading correspondingly to a change of the size and the mass of the configuration under 
consideration. Taking all this into account, we can rewrite (j2"Tj) - (l2~3")) as follows: 



do 

, 2a(n+l)v 

! 1 ~ 6 de 
l + a(n+ 1)0 df 

d^ 



= e 2 <e"^(l + ncr6l)0" + 



1 



1 - 2a(n+l) : 



-e~ H a9 n+1 



2 <r(n + l) 
J ~ l + a(n + l) 



l-2a(n+l)~ 



1 - 2a(n + 1) 



d6 l + cr(n + l)01 ( dv v 



d£ i_ 'Mn+iyv £ £ 



(26) 
(27) 



d£ 



(28) 



These equations are to be solved for given er, n, and b subject to the boundary conditions in the vicinity of the center 
of the configuration £ = 0, 



« 3 r 



(29) 



where cf>o corresponds to the initial value of the scalar field (f>, 6q is normalized to be unity at the center, 9q = 9(0) = 1, 
the parameters 92 , «3 are arbitrary, and the value of the coefficient 02 is defined from Eq. (|28|) as 



h = g b at 



— brt>onn+l 



It is convenient to rewrite the mass and radius of the configuration in terms of the constants K, n, and a (see |21|). 
By eliminating p gc from the expressions for £ and v in (|20p. we obtain 



M(r) = MV (3 - n)/2 u((), M = MV (3 -" )/2 t)fe; 



(30) 
(31) 



where 



i?* = (47T)- 1 / 2 (n + l) 1 /2 G -V2_ ftr n/2 c l-n j 

M* = (4 7 r)- 1 / 2 (n + l) 3 / 2 G- 3 / 2 ^"/ 2 c 3 -". 



(32) 
(33) 



The quantities R* and M* define the scales of the radius and mass. Next, the metric component gn 
can be rewritten via the dimensionless variables (|2"0")) as follows: 



from JT 



—A 



= l-2a(n+l)~. 



(34) 



From Eqs. ([TT]) and (|3"4")l one sees that e u < 1 and e A > 1 along the radius of the configuration, and that e A = 1 at the 
center, while e v takes its minimum value e Ua . Outside the star, the solution approaches the Schwarzschild solution 



= 1 - 2a{n+ 1) 



(35) 



with the total mass M(R) - u(^) [see (p0]l]. 

One may question the assertion that the asymptotic behavior of ([33]) is Schwarzschild since in Sec. IIII Bl we show 
that the scalar field approaches a constant value in the region outside the star. Thus, it would seem that the solution 
should be matched to de Sitter space-time in the external region. However, even though the scalar field approaches 
a constant value in the region outside the star, the energy density of the scalar field goes to zero (see the numerical 
results in Fig. |T|). This comes about since the scalar potential is zero, V(4>) = 0, so that the only contribution to the 
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scalar field energy density comes from the kinetic energy term which involves derivatives of the scalar field. For a 
constant scalar field this kinetic energy term vanishes. 

It will be shown below, in Sec. IIIIB| that the main part of the mass of the configuration under consideration is 
concentrated inside the radius £1. This allows one, by equating expression (|35p to expression (|17[) at the boundary 
of the configuration £ = £1, where 9 = 0, to find the approximate value of the metric function v at the center of 
configuration via 

e- « i-M» + iMft)/a (36) 

[l + a(n + l)] 2 



A. Energies of the system 

In this subsection we review definitions plj ] for the proper energy, the total energy, and the binding energy and 
apply these to the spherical configurations discussed in the previous sections. The purpose is to get some insight into 
the stability of the configurations we will find numerically in the next subsection. In particular, if the binding energy 
of the configuration, as defined below, is negative, then the configuration is naively unstable against dispersal of the 
fluid. A more rigorous, dynamical stability test will be studied in Sec. HVl First, using (pTTj). the total energy E of the 
system, including the internal and gravitational energies, is given by 

E = Mc 2 = M*c 2 a i3 ' n)/2 v{^). (37) 

Second, the proper energy density E is the sum of the rest energy density E 0g of the system, 



E 0g = AI 0g c 2 = 4tt ^ (fPgC 2 + \eT V 2 ) e x ' 2 r 2 dr 



(38) 



and its total microscopic kinetic energy, 



r-R 

£ofc=47r/ fnpe x/2 r 2 dr. (39) 



Using expressions (|T4|) . (fT5j) . and (p4|) for p g , p, and e A , respectively, we get 



£ 0g = Ar c V 3 -")/ 2 Mg (a), (40) 

E ok = A/*cV 3 -")/ 2 Ufc (6), (41) 



where u g and u k are solutions of the equations 



du g _ ^ 2 e- b H" + 1/2 [1 - 2a(n + l)v/tj} (dc^/dQ 2 _ 
d £ [1 - 2<r(n + l)v/£] 1/2 

du k nae-^^e^ 1 



d£ [l-2n{n+l)v/i} 1/2 '' 



u k {0) = 0, (43) 



and Wfc(^i) are the values of these functions at the boundary £ = £i. 

Finally, the binding energy (B.E.) is the difference between the energy of the unbound particles dispersed to infinity 
and the total energy of the bound system, 

B.E. = E Qg - E. (44) 

In terms of u g and v the binding energy can be expressed as 

B.E. = M*cV 3 -")/ 2 [ % (&) - . (45) 

A negative binding energy indicates that the configuration is unstable against the dispersal of the star's matter 
to infinity. Below, we will use this condition to address the question of the stability of the configurations under 
consideration. 
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B. Numerical results 



We solved the system of equations (I26l) - ([28|) numerically using the boundary conditions (|29|) . We started the 
solutions near the origin (i.e. near £ m 0) and solved out to a point £ = £ 1; where the function 8 became zero. Since 
our spherical fluid configuration was taken to be embedded in an external, homogeneously distributed scalar field, we 
also needed to require that at the boundary, £ = £i, the value of the varying scalar field from the inside of the star, 
</>(£), matched the external value </> ex t; i.e. we required 0(£i) — C xt- Under these conditions the numerical calculations 
indicate that on the boundary of the configuration the derivative (dtfi/dt;)^, although small, is not strictly equal to 
zero. It requires matching of the internal solutions with the external ones having a nonzero scalar field energy density 
that vanishes asymptotically. That is why beyond the point £ = £i we continued the numerical solutions with only 
the gravitational and scalar fields while the fluid was set to zero. 

For this purpose, we write down the Einstein equations of Q and (JS|), and the scalar field equation of ([7]), without 
the fluid source, i.e. 8 = 0. This leads to the following system of equations: 



— A 



-A 



12 







+ 1 


(r 2 










+ 1 




4) 




D 2 







AttG 



-e-V 2 , 



AttG 



-e ip 



(46) 
(47) 
(48) 



where D is an integration constant. This system, by analogy with the transformations made above, can be rewritten 
in terms of the dimensionless variables </>(£), an( i ^(0 as follows: 

(49) 
(50) 
(51) 



dv 


ID 2 


-V 


d£ = 




? 


dv 


1 


"1 + 


a(n + l)^e- 


~dl = 


I 


1 - 


- 2a(n + l)v/£ 


r- 


D 2 


e- v 






1 - 


2a(n + l)v/i ' 



where the new dimensionless constant 



D 



AttGD 



a(n + l)c 3 

is introduced. This constant can be defined from (ISTl) as 



D 2 = £ 



l-2a(n+lM^)/^] W 2 } 



Si- 



The system (|49 |) -(|5T |) contains the parameter a as a trace of the influence of the fluid on the external solution. The 
solution is to be sought beginning from the surface of the fluid at £ = £i by using, as the boundary conditions, the 
values of v(£i), <f)(£i) obtained from the solution of equations (|26|) - (f^5|) and from expression (TT71) for the internal 
part of the configuration. This allows one to determine the value of the integration constant v c from (|17[) by requiring 
e v to be equal to unity at infinity, providing asymptotic flatness of the space-time. Thus, the complete solution for 
the configuration under consideration is derived by matching of the internal fluid solutions given by Eqs. (|26l) - (f^H|) 
with the fluid-free, external solutions obtained from the system (|4l^ - (|5"T|) . 

In principle, the total mass of the configuration under consideration is determined as a sum of the internal mass, 
defined by expression (pTTj) . and of the mass associated with the external scalar field. However, in the case of a massless 
scalar field, a more elegant way is to define the total mass via the Komar mass. The latter, in general, is defined 
as 23] 



M. 



K 



r , ; yT ab --g ab T)n a edV, 



where n a is a normal to E and £ is a timelike Killing vector. For the static spherically symmetric configuration being 
considered here, we find 



M K = %\ f(e + 3p)e^ +x ^ 2 r 2 dr. 
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FIG. 1: The distributions of the total energy density in units of p c c 2 from (|25p for the different values of C xt = —2,0,2, from 
top to bottom. At oxt = —2 the central values of T® / p^. are as follows: for n = 1.5, Tq / ' p\ = 8.91; for n = 3.0, T§ / p^. = 10.01 
(not shown in the figure). For all curves 6 = 1. The dashed line corresponds to the case without the chameleon scalar field 
considered in plj. In the n = 1.5, a — 0.1 case of the left figure the energy density smoothly (but rapidly) goes to zero beyond 
the surface of the star as occurs in the right figure. This is not seen explicitly due to the scale used in the left figure. 



Note that here the integration is performed only in the range from to R, where there is a nonzero contribution 
associated with the fluid. Using the dimensionless variables (j2T))) . and also taking into account (|T^|) -([T5 |) . this expression 
can be rewritten in a dimensionless form as 

v K = -^—M K = / ?1 e~ b * [1 + a(n + 3)6} 9"e^ +x ^ 2 ^dt 

Comparing this total mass with the mass of the internal part of the configuration, defined by (|3ip . one can see the 
contribution to the total mass from the external scalar field. The numerical calculations presented below indicate that 
the value of this contribution can be up to 1% in some cases, but generally is much smaller. From this it is clear that 
the main part of the mass is concentrated inside the radius £i (or R). This, in turn, allows one to introduce some 
effective radius of the configuration which is defined as a radius where a certain percentage of mass is concentrated 
(see e.g. Sec. II D in Q). In our case it is natural to choose the effective radius as a radius of the fluid, where more 
than 99% of the total mass of the configuration is concentrated. 

The procedure for finding solutions is as follows: for given values of b, a, and n, we choose the central value of 
the scalar field, <f>o, such that 8 goes to zero and 4>\ ~ (j) cx t at some finite value of £ = £i. As mentioned above, 
£ = £i corresponds to the surface of the star. If one assumes that the value of the scalar field ox t changes during the 
course of the evolution of the Universe, then, in order to study a configuration embedded in such an evolving field, 
it is necessary to specify how the function <j) ex t(t) changes with time. Obviously, this time dependence will be model 
dependent. Here we will use the dependence suggested in [16| . From some general considerations about cosmological 
evolution, the authors of [l6[ showed that the chameleon scalar field may evolve according to the law 

0(t)~ki— — (52) 

where tn is the time when the Universe ends its evolution in the big rip singularity, and the arbitrary factor in front of 
the logarithm we take to be of order unity. In this model the value of the scalar field changes from — oo at the moment 
of the big bang up to +oo at t = i^. It crosses the zero value at t = tji/2. This is called crossing the phantom divide 
since the effective equation crossed from greater than w e ff = — 1 to less than W e ff = — 1. A field with W e ff < — 1 
corresponds to phantom dark energy. This moment corresponds approximately to the present stage of evolution of 
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Relativity parameter log, Q a 



FIG. 2: The dimensionless mass given by Eq. pip versus log 10 a for n = 1.5, 3.0; 6=1; and <j> C xt = —2, 0, 2. The dots represent 
the upper value for stability at log 10 a cr ~ —0.69 for all <^ cx t. 

the Universe. In addition, for the model of [l6|, it is possible to choose parameters so that the coupling function f((f>) 
is positive definite. 

For the same values of the parameters b, a, and n, the central value of the scalar field, <j)Q, will change depending on 
the value of </> ox t. Obviously, this will effect both the radial matter distribution of the star and the star's total mass. 
In Figs. Q][3] the results of the numerical calculations for the three values of (j> cx t — —2, 0, 2 are presented (hereafter we 
choose 6=1). From (|52[) the field values <f> ex t — —2, 0, 2 correspond to the moments of time i_2 ~ 0.12t^, to = 0.5t^, 
and £2 ~ 0.88£t?. The time to corresponds, approximately, to the current stage of the evolution of the Universe, t-2 
to the early Universe, and ti to the future Universe near the big rip singularity. Thus, at different moments of time 
in the evolution of the Universe, the ordinary, gravitating, polytropic matter with fixed parameters er and n creates 
configurations having different distributions of matter along the radius (see Fig. [1} and masses (see Fig. HI). On the 
other hand, these results may be interpreted as an evolution of the physical characteristics of the same star whose 
matter interacts nonminimally with the surrounding chameleon scalar field. 

It is important to note that the cosmological time in (|52p is not the same as the Schwarzschild time of the metric 
([3]). We are assuming that the cosmological evolution is very slow on the time scale for the life cycle of a star. Thus, 
the scalar field is taken to be approximately constant during the life cycle of a given star. In our analysis the only 
dynamical time is the time coordinate from (|3]). The cosmological time coordinate from (|52[) is given to motivate our 
use of different (approximately constant) values of the external scalar field. Of course there could be long-lived stars 
which have a life cycle long enough so that they would need to take into account the effect of the change of the scalar 
field with respect to cosmological time during the life cycle of the star. This is a more complex problem, and we will 
return to it in future work. The general picture that emerges from the /p c c 2 is that the material of the star is more 
concentrated at early cosmological times. This would tend to make these stars hotter and thus have faster life cycles. 
This fits with our assumption that the scalar field is approximately constant during the evolution of these early stars. 
For present times the effect of the scalar field would be minor. This is born out by the fact that the middle solid 
line of Fig. Q] is very close to the dashed line, which is the case without any scalar field. The late time stars have a 
dispersed T® /p c c 2 which would tend to make the stars cooler and thus live longer. It is for these late time stars that 
it might be important to take into account the evolution of the scalar field during the life of the star. 

From Fig. [1] one can see what effect the scalar field has on the configuration. For negative values of cx t [which 
from (|52p correspond to early cosmological times] the effect of the presence of the scalar field is that it tends to 
concentrate the matter of the star with respect to the case without a scalar field. For small values of cx t [which from 
(1521 correspond to the present cosmological time] the scalar field has a minimal effect on the mass-energy distribution 
of the matter. This can be seen by the very small difference between the middle solid line in Fig. [1] and the dashed 
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FIG. 3: The dimensionless binding energy given by Eq. (I45[l versus log 10 a for n = 1.5,3.0; 6 = 1; and </> C xt = —2,0,2. All 
the binding energies eventually become negative, but for n = 1.5 there is first a peak in binding energy after which instability 
occurs. For n = 3.0 the binding energy is always negative. The dots represent the upper value for stability at log 10 a C r — —0.69 
for all 0ext- 

line which corresponds to no scalar field at all. Note, even though <fi ext = for the middle solid line in Fig. [TJ there 
is a small difference between this situation and having no scalar field, since even though </> ext = there is still some 
nonzero distribution of the scalar field inside the star which leads to this small difference. Finally, for positive values 
of (frext [which from (|52|) correspond to late cosmological times] the scalar field tends to diffuse the matter of the star. 
Thus, the general picture that emerges is that the effect of the scalar field at early times is to concentrate the matter 
of the star, which would have the effect of making these stars hotter than without the effect of the scalar field; at the 
present times the effect of the scalar field would be minimal; at late times the effect would be to disperse the matter 
of the star, which would have the effect of making these stars cooler than without the effect of the scalar field. 

Figure [2] shows the dependence of the dimensionless mass M/M* as a function of log 10 a for two values of the 
polytropic index n — 1.5 and n = 3.0, and the values of the scalar field at the surface of the star, 4>i ~ cxt = —2, 0, 2. 
Here we have restricted ourselves to just these two values of n since, as in the case of the polytropic sphere without 
a scalar field considered in plj . they provide qualitatively different behaviors of the mass as a function of a. In the 
case n — 1.5, as opposed to the case n = 3.0, the mass as a function of a rises to a maximum and then undergoes 
damped oscillations about an asymptotic value. 

Similar differences in the behavior of the configurations at n = 1.5 and n = 3.0 can be seen from Fig. |3l where 
the dependence of the dimensionless binding energy B.E. /M*c 2 as a function of log 10 a is shown. For n — 1.5 there 
is some value <7/ p of er at which the binding energy has its first peak. Configurations having a < <7f p have positive 
binding energy, and for a > <7f p the binding energy decreases and eventually becomes negative. In turn, for n = 3.0 
the binding energy is always negative. 

As mentioned above, the negativeness of the binding energy implies that the configurations will be unstable against 
the dispersal of the fluid of the star to infinity. However, even having a positive binding energy does not guarantee 
the stability of a configuration. For example, with n = 1.5, and for the configurations considered in [2l|, the binding 
energy remains positive for a > af p , but the configurations are unstable - it is energetically favorable to make a 
transition to a state with the same polytropic index n, but with a smaller parameter a. Thus, the first mass peak for 
the case n = 1.5 from [2l| corresponds to the point dividing stable and unstable configurations, and in this case the 
value afp corresponds, in Tooper's terminology, to the critical value a cr . For the case of n = 3.0 there is no peak in 
either Tooper's model or in the problem being considered here: all configurations have negative binding energy and 
are unstable for all a and cx t- 

The question arises as to whether a f p coincides with a cr in the presence of the chameleon scalar field for the case 
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n = 1.5. In |2Jj this question was studied for a boson configuration supported by the complex scalar field using the 
fractional anisotropy of the system. The fractional anisotropy (fa) is defined by 

fa = (p r - pt)/p r , (53) 

where the radial pressure p r = —T\ and the tangential pressure pt = — Tf = — Tg are given by the corresponding 
energy-momentum tensor. Using the variational principle suggested by Chandrasekhar [25[, Gleiser [24j applied the 
method of trial functions to estimate the dynamical stability of this system. As a result, contrary to the results for 
fluid stars, Gleiser found that the instabilities did not occur for central densities which were equal to or larger than 
the critical density but for configurations well to the right of the maximum mass. However, further studies of such a 
system showed that the trial functions used in [24| were not the most general ones, and that the instability occurred 
at the maximum mass, as in the case of ordinary neutron stars (see the erratum |24j|). 

However, in our system, in addition to the fractional anisotropy, we also have nonminimal coupling between ordinary 
matter and the scalar field. Because of this difference, in the next section, we will consider in detail the question of 
dynamical stability of the configurations studied here in order to clarify the possible influence of such nonminimal 
coupling on the stability of the system. In doing so, the critical value a cr will be calculated by studying the dynamical 
stability of the configuration using the approach suggested by Chandrasekhar [25j]. We will not use the method of 
trial functions used in [24j . but we will instead numerically integrate Chandrasekhar's pulsation equations, to trace 
the behavior of the lowest frequency modes and thus determine a cr . 



IV. DYNAMICAL STABILITY 



A. General equations 

Consider that the equilibrium configurations described in the previous section are perturbed in such a way that 
spherical symmetry is maintained. In obtaining the equations for the perturbations, we will neglect all quantities 
which are of the second and higher order. The components of the four-velocity are given by 25] 

—Uo/2 i/n/2 1 -v a /2 Ao-i/n/2 



with three- velocity 



dr 



The index in the metric functions indicates the static, zeroth order solution of the Einstein equations. The compo- 
nents of the energy-momentum tensor ^ are then given by 



n = f* + \e~ V + \e~ V 2 + V(<p), 

Tl = -fp V V - \e~ V 2 + 

t q = f( £ +p)u u 1 + <9 <ys<9V = /(eo +Po)v - e~ X ipip' 7 

2 rji3 



-fp-\e-^ 2 + \e-\' 2 + V^). 



Now we consider perturbations of the static solutions of the form 

y = yo + y P , 



(54) 

(55) 
(56) 
(57) 

(58) 



where the index refers to the static solutions, the index p indicates the perturbation, and y denotes one of the 
functions X,u,e,p or (p. Substituting these perturbed expressions into Eqs. Q and ([5]), we find 



-Ao 



d 



[rX' p + \ p (l-r\' )]^-[re-""X p ] 



8irG 



dr 

e- Xo [rv' p - X v (1 + rv' )] = ^r 2 f oPp + f pPQ + e~ A Vo W P ~ \<Ja>*) 



- A o,„' 









(4 




+ v p 


-v P 







(59) 
(60) 
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where 

fo = f(<Po), fp = <P P 
Next, from Eq. (|6|) we have 



dp J ip=ip 



, V Q = V(ipo), V p = <p p 



dV 

dip J ip=ip 



- Xo ^ = ^f-[f (£ +p )v-e->° i p p <p' ] 



Using the components (|54l) - (|57[) in ([8]l. we can rewrite this equation in terms of the perturbations as 
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[/0O0 + Po)v - e °ip' <Pp] - q- (foPp + fpPo) 



<9r 



fo(£p+Pp) + f P (£o+Po) v' Q - -fo{e +PoW p 



Now it is useful to introduce a "Lagrangian displacement" ( with respect to x°, 



Equation (jBTjl can be directly integrated to give 



-^eW[f (e a +p a )C 



Inserting this expression into Eq. (I59[) gives 



1 



Jo 



y + — — 

P r 2 dr 



r 2 (/o(£o+Po)C-e- A Vo^) 



The expression for vL follows from Eq. (pD|) . 

v'={1 + n/ ) + — j-» 



foPp + fpPo + e °v'o[v'p- -^f'aK) - y p 



The perturbed scalar field equation can be found from (J7J, and is 



^ - e V P 



dV\ 
dip J, 



Pa 



df_ 

dip 



P P 



df 



+ A p 



where 



dip 



d(f,V) 

dip 



dip J 
d(f,V) 



dV 



Po 



tp=ip 



dp 



dtp J 

d 2 (f,V) 



dip J . 



Vp 



dip 1 



<p=ip 



(61) 



(62) 



(63) 



(64) 



(65) 



(66) 



Thus, we have two linear second-order partial differential equations (|62[) and (|66p for the displacement £ and the 
perturbation of the scalar field ip v . To write out these equations, one needs to use expressions (|63[) - (|65[) . 
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B. 



Case of f — e 



In this section we apply the equations obtained in the previous subsection to study the stability of the spherical 
solution considered in Sec. IIIII In this case V = 0, and the equation of state of matter is taken from expressions 
(I13I) - (I15I) . From these equations one finds that the perturbed components of the pressure p p and the energy density 
s ji are 



Pp = K(n + l) P ]t 1/n e^9 p , e p = n Pgc c 2 
In turn, the static components are 



1 



+ a{n + 1) 



Po = Kp)+ l/n 0% + \ £ = p gc c 2 (1 + an9 ) 0J. 
Next, using expressions ([57]) and (j6"5)l and the dimensionless variables from (f2"0")) we write Eqs. (p5|) -(|55 |) as 

A p = -2a(n + l)e A °£ {f [1 + a(n + l)9 ] 9^ ~ e~ Ao 0o<M , 

x { [/„(! + anOo) + /o + «r(n + l) 2 ] tftfl) + [1 + a(n + 1)0 ] (fo^ + ~M + 



(67) 

(68) 
(69) 



,-Ao 



(fo + %fo - K^'o)^ + \^\ }> 



v' P = (1 + y + Mn + l)e x ^[a[f (n + 1) 



'o w p 



fp 9 o 



n+1 



,-An J.' 



(70) 
(71) 



where / = e" 6 ^, / p 



and we have introduced the dimensionless displacement ip — AC- Equations 



([6^| - ([7T|) will be used to turn the general equations (j6"2")l and (|66p toward the specific nonminimal coupling under 
consideration here, namely, / = e~ b ^. 

To proceed with the stability analysis we assume harmonic perturbations having the following time dependence, 



(72) 



The function y p (£) depends only on the space coordinate £. For convenience, we hereafter drop the bar. Using (|72j) . 
the general equation (|62[) takes the form (we have introduced a new dimensionless frequency ui = lo/A and we again 
subsequently drop the bar for simplicity) 



-An Aj 



u 0Vp} 



,-Ao 



[fo [1 + a(n + 1)6 ] 9^ 
| (fo - K<P'o) (fp ~ t^o a p) + 0p' - 5 ^" Ap + ^ A p) 



/o(£ + <r(n + 1) 2 )^ P + / P (l + a(n + l)^o)^] ^ - 5/0(1 + + l)0 )e^' p 



(73) 



-Ao 



-00 + O U> p - -0q A P v 



0. 



Next, the perturbation equation for the scalar field (|66]l is given by 



1 



(74) 



= -(7e Ao {6 2 /oCVp " bfo [(n + 1)9 P + 9 \ P ] 0%) 



Thus, we have two linear second-order ordinary differential equations ([73"]) and (j?4"|) for the dimensionless displacement 
ij) and the perturbations of the scalar field <f> p . Equation (|T3[) is an analog of Chandrasekhar's "pulsation equation" 
25]. Solutions of this equation must satisfy the boundary conditions ip(0) = and p p = at £ = £i. On the other 
hand, since the total expression for the perturbed radial pressure Sp r = —8Tl [with T-j 1 taken from (|55p] contains 
not only the fluid component p p but also perturbations of the scalar field, it is necessary to require that the total 
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FIG. 4: The perturbation in the metric function A p is shown as a function of the radial coordinate £ for the first three modes 
for the critical value log 10 a cr ~ —0.69. The curves are plotted for 6 = 1 and </> cxt = 0. 

perturbation 5p r be equal to zero as well. Formally, the scalar field energy goes to zero at infinity, i.e. as £ —> oo. 
Thus, the boundary conditions for the system of coupled equations (|73|) and (|74p are 

ip(Q) = <f> p (0) = 0, S Pr -> as £ ^ oo. (75) 

The system ((75)1 and (f74)) . together with the boundary conditions ((75)) . defines a characteristic value problem for cj 2 . 
The question of stability is thus reduced to a study of the possible values of uj 2 . If any of the values of uj 2 are found to 
be negative, then the perturbations will grow and the configuration is unstable against radial oscillations. Since the 
eigenvalues form a discrete sequence, our purpose is to find the lowest eigenvalue ui 2 , and if uj 2 , > the configuration 
is stable. 

Of course, the value of ujq depends on the values of the parameters n, a, b, and the value of the background scalar 
field 4> cx t- Here we find the values of a for which the equilibrium configurations are stable at fixed n, b, and </> cx t- This 
will allow us to determine the boundary between stable and unstable configurations. This boundary is defined by the 
point where uj 2 ., as function of a, crosses through zero. We denote the corresponding critical values of a as a cr . 

To find cr cr , we numerically solved Eqs. ((73)) and (f74| to determine the characteristic values for the frequency uj 2 . 
In the absence of the scalar field, both the energy approach to stability and the dynamical stability analysis give the 
same values of a cr [21[. Further, without a scalar field, a cr is given by the first peak in both the mass as a function of 
a and the binding energy as a function of a distributions. However, in the presence of a scalar field the point dividing 
stable and unstable configurations may shift. Therefore, in finding er cr , we proceed as follows: we choose the initial 
value of a corresponding to the first peak in the binding energy (or in the mass) functions. For this value, we find the 
lowest 0Jq. We always get positive values of uj 2 , at the first peak, which means that these configurations are stable. 
Next, shifting to the right from the first peak, we look for the value of a at which uJq = 0, which gives a cr . 

We applied the above procedure for finding a cr for the configurations with the fixed values of n, b, and 4> ex t used in 
Sec. lIIIBl The values of a cr are shown in Figs. [2] and [3J by the bold dots. The dots mark the dividing point between 
stable configurations (to the left of the dots) and unstable ones (to the right of the dots). As expected, for the case 
n = 3.0, when the binding energy is always negative, there are only negative values of ujq; i.e. the configurations are 
unstable at any a. Thus, in Figs. [2]and|3j only the results for the case of n = 1.5 are given. In this case the numerical 
calculations give log 10 a cr ~ —0.69 for all values cxt = —2, 0, 2. In turn, the fractional anisotropy (|53p is given by 

fa= [l-2<x(rc+l)v /flfl, 2 

vM£ +1 + 1/2 [1 - 2a(n + l)v /£] <t>'$ ' 

It can be found from this expression that at the center of the configurations fa(0) = 0, and at the boundary, £ = £i, 
fa(£i) is always equal to 2. 

The eigenfunctions corresponding to the lowest Wq have no nodes, and the eigenfunctions corresponding to the nth 
mode have n nodes. As an example, Fig. U shows the plots of the eigenfunction A p from (|69[) for the above value 
log 10 a cr ~ —0.69. Using Eqs. ((73")) and ((74")) and the boundary conditions ((75)) . we plot in this figure the first three 
modes. For the zero mode ujq = 0, for the first mode uj\ ~ 0.18, and for the second mode uj\ ~ 0.22. 
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V. CONCLUSION 

Based on the assumption that time evolving chameleon scalar fields may exist in the Universe, this paper studies the 
possible influence that such fields may have on the inner structure and stability of compact gravitating configurations 
consisting of ordinary matter and embedded in an external, homogeneous chameleon scalar field. The ordinary matter 
is taken to have a polytropic equation of state which is parametrically given by expressions (p~3|) - (p~5|) . We choose 
a specific form for the nonminimal coupling between the ordinary matter and the scalar field given by (1241) . This 
form was used in [TH ] within the framework of chameleon cosmology. Here we study the consequences of the presence 
of a chameleon scalar field on the inner structure of compact astrophysical objects (i.e. "stars") rather than the 
cosmological consequences of chameleon fields. In Ref. [TtJ this type of gravitating configuration of a polytropic fluid 
nonminimally coupled to a scalar field was called a "chameleon star." 

One of the main conclusions of the present work is that the characteristics of the chameleon stars - their masses, 
sizes, matter distribution - depend strongly on the parameters of the surrounding scalar field <f>. Assuming that this 
background scalar field changes with cosmological times, we predict that the characteristics of the stars embedded in 
the scalar field should also change. For example, we considered three different values of the background scalar field, 
0cxt = —2, 0, 2, which could represent the time evolved value of the field at different times according to the law (|52p . 
These three values could correspond to three different periods of time in the evolution of the Universe - from the 
relatively early Universe to the future Universe, for example, already near the big rip singularity. The results of the 
numerical calculations presented in Figs. [T][3] (which cover the values of the polytropic index, n = 1.5 and n = 3.0) 
indicate that the characteristics of the configurations have a strong dependence on the value of the scalar field <fr cx t- 
In particular, from Fig. [T]one finds that, at the center of the configurations considered, there can be either a larger or 
smaller concentration of the matter density compared to the case of a star without a scalar field. This would have an 
impact on the rate of fusion reactions inside the stars, and thus on their lifetimes so that the life span of the star could 
be altered by the properties of the background chameleon scalar field surrounding the star. Some other speculations 
about possible physical applications of the chameleon star model can be found in Ref. [I?} ■ 

Another important consequence of the presence of the chameleon scalar field is its influence on the stability of the 
configurations studied here. From previous investigations of the stability of polytropic gravitating fluids it is known 
that for equilibrium configurations instabilities occur once one moves to the right of the first peak in the binding 
ener gy ( or the first peak in the mass) [2lj . A similar situation also occurs for purely boson stars (see the erratum from 
Ref. [24j|). However, the dynamical stability analysis performed in Sec. II VI shows that the configurations considered 
here have their stable points shifted to the right of the first peak in the binding energy. These shifts are shown by 
the bold dots in Figs. [5] and [3] Thus, in contrast to the usual polytropic and boson stars, the expansion of the region 
of stability is in the direction where a increases. It is obvious that this new effect is caused by the presence in the 
system of the nonminimal coupling between the ordinary matter and the scalar field. 

The study of mixed configurations with scalar fields and ordinary matter sources is of interest since various types 
of scalar fields are believed to exist in the Universe. Most often, such scalar fields play a cosmological role - as a 
mechanism to drive inflation (the inflaton field), as dark matter, or as dark energy. Here we consider the possibility 
that one type of scalar field (a chameleon scalar field) might also have some influence on astrophysical objects such 
as stars. 
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